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Abstract. In this talk I review various developments in the past year concerning quantum 
vacuum energy, the Casimir effect. In particular, there has been continuing controversy 
surrounding the temperature correction to the Lifshitz formula for the Casimir force between 
real materials, be they metals or semiconductors. Consensus has emerged as to how Casimir 
energy accelerates in a weak gravitational field; quantum vacuum energy, including the divergent 
parts which renormalize the masses of the Casimir plates, accelerates indeed according to the 
equivalence principle. Significant development has been forthcoming in applying the multiple 
scattering formalism to describe the interaction between nontrivial objects. In weak coupling, 
closed-form expressions for the Casimir force between the bodies, which for example reveal 
significant discrepancies from the naive proximity force approximation, can be achieved in many 



1. Introduction 

There has been burgeoning interest in the Casimir effect [1], the quantum force between 
macroscopic bodies due to zero-point fluctuations, as witnessed by the multitude of international 
conferences and workshops devoted to the subject, including two in Brazil back-to-back. This 
interest reflects a great deal of progress in the subject both from the experimental and theoretical 
sides. I will talk about three topics my collaborators and I have been involved with during 
the past several months, which mirror these developments: the controversy concerning the 
temperature dependence of real materials, such as metals and semiconductors; the question of 
how the Casimir energy associated with parallel Casimir plates is accelerated by gravity; and the 
achievement of exact closed-form results for the interaction of bodies subject to weak coupling. 



2. Temperature Controversy 

2.1. Introduction 

There is a continuing controversy surrounding the question of how to incorporate thermal 
corrections into the Casimir force between real metal plates [2,3]. The procedure to do this 
seems straightforward. The thermal Green's functions must be periodic in imaginary time, with 
period /3 = 1/T. This implies a Fourier series decomposition, rather than a Fourier transform, 
where we have a sum over Matsubara frequencies Cm = 2irmT. 

1 The work reported here was done in collaboration with I Brevik, S Adn0y Ellingsen and J S H0ye (section [2}; 
S Fulling and A Romeo (section [3]) ; I Cavero-Pelaez, P Parashar, K V Shajesh and J Wagner (sections [3] and UJ). 



Thus, the Casimir (Lifshitz) [4] pressure between parallel semi-infinite nonmetallic plates 
separated by a distance a is given by 

t=TE,TMm=0 J <> m y J 

where the prime means that the m = term is counted with half weight, in terms of the reflection 
coefficients 
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Let us rewrite the Lifshitz formula at finite temperature in the form 
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where the second equality uses the Euler-Maclaurin sum formula, in terms of 
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where 

4n = (rT 2 e 2Ka ~ 1, (5) 

and k 2 = kj_ + C 2 - 

Evidently, for the Drude model, where the permittivity has the form 



or more generally, whenever 
the TE zero mode vanishes, 
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This leads to a physical discontinuity in the TE mode, between the m = and m = 1 Matsubara 
frequencies. Naively, this implies an additional linear term in the pressure at low temperatures: 

P T „ pT=0 + 03) ^ <<; L (g) 
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Exclusion of the TE zero mode will also reduce the linear temperature dependence expected at 
high temperatures, 

pT ~"!^l T > aT»l, (10) 
otto 13 



one-half the usual ideal metal result. This anomaly was first pointed out by Bostrom and 
Sernelius [5]. See also [6]. 

Several other recent papers also lend support to our point of view. Jancovici and Samaj [7,8] 
and Buenzli and Martin [9, 10] have examined the Casimir force between ideal-conductor walls 
with emphasis on the high-temperature limit. Not surprisingly, ideal inert boundary conditions 
are shown to be inadequate, and fluctuations within the walls, modeled by the classical Debye- 
Hiickel theory, determine the high temperature behavior. The linear in temperature behavior of 
the Casimir force is found to be reduced by a factor of two from the behavior predicted by an 
ideal metal, just as seen above in (|lUp . This is precisely the signal of the omission of the m = 
TE mode. Further support for our conclusions can be found in the work of Sernelius [11], who 
calculated the van der Waals-Casimir force between gold plates using the Lindhard or random 
phase approximation dielectric function. The central theme of his work is to describe the thermal 
Casimir effect in terms of spatial dispersion. He also finds that for large separations the force 
is one-half that of the ideal metal. Sernelius shows that, for arbitrary separation between the 
plates, the spatial-dispersion results nearly exactly coincide with the local dissipation-based 
results. 

Thus, it is very hard to see how the corresponding modification of the low-temperature 
behavior can be avoided. However, because a linear T dependence in the pressure implies a linear 
T dependence in the free energy, this seems to imply a violation of the Nernst heat theorem, 
the third law of thermodynamics, because then a naive calculation implies S{T = 0) ^ 0. 
However, we [12] and others have shown that for real metals, where 7 has a residual value at 
zero temperature, the free energy per area has a vanishing slope at the origin. Indeed, in the 
Drude model, the free energy per area has the behavior [13, 14] 

LO 2 

F = F + T 2 ^-(21n2-l), (11) 

for sufficiently low temperatures. The T 5//2 and T 3 corrections to this were calculated last year 
in [12,13] This is illustrated in figure [U taken from [12], which shows the numerical evaluation 
of the free energy for gold. 

2.2. Experimental constraints 

The difficulty is that, experimentally, it is not easy to perform Casimir force measurements 
at other than room temperature, so current constraints on the theory all come from room 
temperature experiments. Then all one can do is compare the theory at room temperature 
with the experimental results, which must be corrected for a variety of effects, such as surface 
roughness, finite conductivity, and patch potentials. There are no direct experiments which 
observe T dependence in metals! (The temperature dependence of the Casimir-Polder force 
between a Bose-Einstein condensate and a dielectric substrate was measured in 2006 [15].) 

Although the earlier experiment by Lamoreaux [16] was carried out at such large distances 
(~ 1 jum) where the temperature correction proposed here would be large (~ 15%) [17]. there 
are reasons to think that the uncertainties due to such things as surface roughness and patch 
potentials were so large as to afford only limited accuracy in his experiment. There has been a 
great number of high quality experiments during the last decade, revealing a number of facets 
of the quantum vacuum energy [18-27]. As for the temperature effect, the most stringent test is 
given by the precision experiments of Decca et al [28, 29] , where the agreement of the data with 
the zero-temperature theory seems to rule out a large thermal correction such as Bostrom and 
Sernelius [5] and we have proposed. It seems to us that there may be a number of effects that 
have not been properly accounted for before one can conclude that no large temperature effect 
is present, such as the shielding effect recently proposed by Pitaevskii [30]. 
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Figure 1. Numerical evaluation of the free energy between two gold halfspaces as a function of 
temperature. The inset gives details for low T. 



2.3. Anomaly for semiconductors 

Recently Klimchitskaya, Geyer and Mostepanenko [31-34] have observed there is a similar 
anomaly, now affecting the TM reflection coefficient, for dielectrics and semiconductors. Here 
is a simple way to understand that argument. Suppose we model a dielectric with some small 
conductivity by the permittivity function 
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(The only essential point is that as £ — > 0, e — > e if a = 0, otherwise e — * oo.) 

For a dielectric the TE reflection coefficient is continuous, but if there is vanishingly small 
(but not zero) conductivity the TM coefficient is discontinuous: 
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This implies a linear term in the pressure 
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where the poly logarithmic function is 
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Note that the linear term vanishes for e — » oo. The free energy term F is obtained from this by 
multiplying by a/2. Thus at zero temperature, the entropy is nonzero, 
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Moreover, besides this thermodynamical contradiction, recent experiments by Mohideen's group 
[27, 35-37] show that apparently the conductivity effect must be excluded for semiconductors 
with very low carrier concentrations, but included for semiconductors with high carrier 
concentrations. This conundrum has not yet been resolved, although there has been some 
progress in understanding what happens when the conductivity remains small but nonzero 
at zero temperature [38] — see also Simen Adn0y Ellingsen's talk in this meeting [39]. We 
should also mention our joint work on temperature effects and anomalies [40] which, while not 
addressing the semiconductor anomaly, attempts to survey the underpinnings of thermal effects 
in the Casimir effect. 



3. How Does Casimir Energy Fall? 

In a series of papers last year [41-43] we showed that the total Casimir energy of an apparatus 
consisting of parallel plates, including the divergent parts of the energy, which renormalize the 
masses of the plates, possesses the gravitational mass demanded by the equivalence principle. 

This result puts the lie to the naive presumption that zero-point energy is not observable. 
On the other hand, because of the severe divergence structure of the theory, controversy has 
surrounded it from the beginning. Sharp boundaries give rise to divergences in the local energy 
density near the surface, which may make it impossible to extract meaningful self-energies of 
single objects, such as the perfectly conducting sphere considered by Boyer [44]. These objections 
have recently been most forcefully presented by Graham, Jaffe, et al [45], but they date back 
to Deutsch and Candelas [46]. In fact, it now appears that these surface divergences can be 
dealt with successfully in a process of renormalization, and that finite self-energies in the sense 
of Boyer may be extracted [47, 48] . 

Gravity couples to the local energy-momentum tensor, and such surface divergences promise 
serious difficulties. We first ask how does the completely finite Casimir interaction energy of 
a pair of parallel conducting plates, as well as the divergent self-energies of non- ideal plates, 
couple to gravity? Disparate answers had been given in the past [49,50]; thus, the question, and 
its answer, turn out to be surprisingly less straightforward than the reader might suspect! 

We will not here address what gravitational field results from the (in general divergent) 
Casimir energy. For a beginning of the renormalization of Einstein's equations resulting from 
singular Casimir surface energy densities see [51]. 

Brown and Maclay [52] showed that, for parallel perfectly conducting plates separated by a 
distance a in the z-direction, the electromagnetic stress tensor acquires the vacuum expectation 
value between the plates 

( T ^) = ^diag(l,-l,-l,3), £ c = "^3^. (17) 

Outside the plates the value of (T' 11 ') = 0. 

Because there are some subtleties here, let us review the derivation of (|17j) for the case of 
a conformally coupled scalar (the electromagnetic case differs by a factor of two). Actually, 
the result between the plates, < z < a, is given in great detail in [53] (£ is the conformal 
parameter) 

(TH = (n + u)diag(l, -1, -1, 3) + (1 - 6£)s(s)diag(l, -1, -1, 0), (18) 



where 

I f°° 7T 2 
Un = - / dKK 3 , U= 7, (19) 
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Note that no is a divergent constant, independent of a, and is present (as we shall see) both 
inside and outside the plates, so it does not contribute to any observable force or energy (the 
force on the plates is given by the discontinuity of {T zz }), and so may be simply disregarded 
(as long as we are not concerned with dark energy). Similarly, in terms of the Hurwitz zeta 
function, 

1 

which exhibits the universal surface divergence near the plates, 



s ( z ) = [C( 4 » *A0 + C(4, 1 - z/a)}, (20) 



S(Z) ~-16^' ^° + ' (21) 

is also unobservable (if we disregard gravity) because it does not contribute to the force on the 
plates, nor does it contribute to the total energy, since the integral over s(z) between the plates 
is independent of the plate separation. Of course, the best way to eliminate that term is to 
choose the conformal value £ = 1/6. 

Since the exterior calculation does not appear to be referred to in [53], let us sketch the 
calculation here: Consider parallel Dirichlet plates at z = and z = a. The reduced Green's 
function satisfies 

-£ s + iAg(z,z')=S(z-z'), (22) 
where k 2 = k 2 — uj 2 = k 2 + ( 2 . The solution for z, z' < is 

g (z, z') = --e KZ< sinhKz>. (23) 

K 

where z <t> means the lesser or greater of z and z'. It is very straightforward to calculate the 
one-loop expectation value of the stress tensor from 



(TH = (d»d' v - ^ v d x d'^j ~G(x,x') 



i{d"d v -g^d 2 )^G{x,x), (24) 



which follows from 

(<j>(x)<j>(x')) = jG(x,x'). (25) 
After integrating over uj = \Q and k, we find the result (z < 0) 

(T^ v ) = n diag(l, —1, —1, 3) - ^ ~ 6 ^ diag(l, -1, -1,0). (26) 

This is exactly as expected. The no term is the same as inside the box, so is just the vacuum 
value, and the second term is the universal surface divergence seen in (I2ip (independent of plate 
separation), which can be eliminated by setting £ = 1/6. 

Thus, we conclude that the physical stress tensor VEV is just that found by Brown and 
Maclay [52]: 

(T^) = itdiag(l,-l,-l,3)0(z)0(a-z). (27) 



in terms of the usual step function. 
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Figure 2. Relation between two Cartesian coordinate frames: One attached to the earth 
(x,y,z), where — z is the direction of gravity, and one attached to the parallel-plate Casimir 
apparatus (Cj^jX); where £ is in the direction normal to the plates. The parallel plates are 
indicated by the heavy lines parallel to the r\ axis. The x = x ax is is perpendicular to the page. 



3.1. Variational principle 

Now we turn to the question of the gravitational interaction of this Casimir apparatus. It seems 
this question can be most simply addressed through use of the gravitational definition of the 
energy-momentum tensor, 

5W m ^- l -j (dx)V=gSg^T tlv = IJ (dx)v^^T^. (28) 

For a weak field, 

9ij,u = V^f + 2h fMU , (29) 

so if we think of turning on the gravitational field as the perturbation, we can ignore \J—g. The 
gravitational energy, for a static situation, is therefore given by [SW = — J dt 5E) 



E g = - J (cbOV?^. (30) 

Here we use the gravity- free electromagnetic Casimir stress tensor (|27|) . with u now replaced by 
£ c /a for the electromagnetic situation. 

The Fermi metric describes an inertial coordinate system, as we will discuss further below in 
section 13.21 

froo = ~gz, h 0i = hij = 0. (31) 

Let us consider a Casimir apparatus of parallel plates separated by a distance a, with transverse 
dimensions I>a, Let the apparatus be oriented at an angle a with respect to the direction of 
gravity, as shown in figured! The Cartesian coordinate system attached to the earth is (x, y, z), 
where z is the direction of — g. The Cartesian coordinates associated with the Casimir apparatus 
are (£, r/, x), where £ is normal to the plates, and rj and x are parallel to the plates. 
The relation between the two sets of coordinates is 

z = Ccosa + 77 sin a, (32a) 
y = 7/ cos a — Csina, (326) 
x = x- (32c) 

Let the center of the apparatus be located at 

(C = Co,?7 = 0,x = 0). (33) 



Now we calculate the gravitational energy from (|30f) . 

[■ g cL/2 Ko+a/2 

E g = (dx)gz(T 00 ) = —gL drj d((( cos a + rj sin a) = g£ c L 2 ( cos a + K, (34) 

J a J-L/2 JCo-a/2 

where K is a constant, independent of the center of the apparatus Co> z o = Co cos a. Thus, the 
gravitational force per area A = I? on the apparatus is independent of orientation 

a small upward push. Here e = \ga is a measure of the gravitational force relative to the Casimir 
force J- c = — 7r 2 /(240a 4 ). Note that on the earth's surface, the dimensionless number e is very 
small. For a plate separation of 1/um, 

6 = § ~ 10 ~ 22 ' ( 36 ) 

so the considerations here would appear to be only of theoretical interest. The effect is far 
smaller than the Casimir forces between the plates. 

It is a bit simpler to use the energy formula ([30]) to calculate the force by considering the 
variation in the gravitational energy directly, as we illustrate by considering a mass point at the 
origin: 

= m5(r)5^5 u0 . (37) 

If we displace the particle rigidly upward by an amount Szq, the change in the metric is 
^oo = —g5zQ. This implies a change in the energy, exactly as expected: 

5E g = —m (—g5zo) = mg5zQ. (38) 

Now we repeat this calculation for the Casimir apparatus. The gravitational force per area 
on the rigid apparatus is 

F = -f^ = -g£ c , (39) 
oz 

again the same result found in (|35|) . which is 1/4 that found by Bimonte et al [50]. It does, 
however, agree with one of the results found in the earlier paper by the same collaboration [49], 
who now completely agree with our calculation [54]. Our result is further completely consistent 
with the principle of equivalence, and with one result of Jaekel and Reynaud [55]. 
In electrodynamics, one defines the field (the 4- vector potential) by 

5W = J {dx)5J^A^. (40) 

From this, one can compute the force on a charge by considering the displacement of a static 
charge, 

JA» = ert(r), (41) 

that is, 

SJi* = eSt* (Sr • V)5(r). (42) 
Inserting this in the above and integrating by parts we get 

5W = J dt J (dr)e (-VA ) ■ St6(t). (43) 



Because the electric field is E l = F° l = d°A l — d i A°, and a total derivative on time is irrelevant, 
we deduce the force on the charge to have the expected form: 

- = F = eE. (44) 
Or 

We therefore should be able to proceed similarly with vacuum energy, starting from the 
definition of the gravitational field 



SW = J{dx)5T» v h^. (45) 
Again, check this for the force on a mass point, described by ([57]) . so 

5E g = - j (dr)m (-<5r • V) 5(r)h 00 = -m5r ■ V/i 00 . (46) 
since h°° = -gz, we conclude 

5E , , 

- — = F = -mgz. (47) 

For the same constant field the force on a Casimir apparatus is obtained from the change in 
the energy density 

(TOO) = ^e{a/2 - C + (o)9(( - Co + a/2), (48) 
a 

that is, recalling that z$ = Co cos a, we obtain the change of the energy density under a rigid 
displacement 8zq: 

S(T°°) = ^5z — [S(( - Co - a/2) - 6(( - Co + a/2)} , (49) 
a cos a 

which yields the gravitational force, since h 00 = —g((cosa + ^sina), 



jF = SE 9 = = £ c 1 ^oo 
A5zq A a cos a 



C=Co+a/2 

= -g£ c . (50) 

C=Co-a/2 



3.2. Inertial coordinates 

However, one might think that the above metric, while sufficing for massive Newtonian objects, 
might seem inappropriate for photons, near the surface of the earth. Rather, shouldn't we use 
the perturbation of the Schwarzchild metric, for weak fields (GM/r <C 1) either in its original 
or isotropic form? Doing so would give rise to a different force on the Casimir apparatus. The 
reason we get different answers in different coordinate systems reflects the fact that our starting 
point is not gauge invariant. Under a coordinate redefinition, which for weak fields is a gauge 
transformation of h^ v [56] 

hfa, ->• hfa, + d^ v + d v ^, (51) 

where C^ is a vector field, the interaction W is invariant only if the stress tensor is conserved, 
dpT^ = 0. Otherwise, there is a change in the action, 



AW = -2 (dx)C^T^. 



(52) 



Now in our case (where we make the finite size of the plate explicit, but ignore edge effects 
because L»a) 



TV = ^diag(l, -1, -1, 3)6 (c-Co + |)^(|-C + Co) 
x9( V + L/2)8(L/2 - n)6( X + L/2)9(L/2 - X ). 



(53) 



Thus the nonzero components of dJTV are 



[5(C - Co + a/2) - S(( - Co - a/2)] 6. . . 
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-^[5( x + L/2)-5( X -L/2)]6..., 



(54 a) 




(546) 



(54c) 



where 6 . . . refer to the remaining step functions. Therefore, the change in the energy obtained 
from AW is 



Since we have demonstrated that the gravitational force on a Casimir apparatus is not a 
gauge- invariant concept, we must ask if there is any way to extract a physically meaningful 
result. There seem to be two possible ways to proceed. Either we add another interaction, say a 
fluid exerting a pressure on the plates, resulting in a total stress tensor that is conserved, or we 
find a physical basis for believing that one coordinate system is more realistic than another. The 
former procedure is undoubtedly more physical, but will yield model-dependent results. The 
latter apparently has a natural solution. 

A Fermi coordinate system is the general relativistic generalization of an inertial coordinate 
frame. Such a system has been given in [57] for a resting observer in the field of a static mass 
distribution. For our case, where the gravitational potential is <j) = gz (apart from an irrelevant 
additive constant), we recover the Fermi coordinate metric (]3ip for a gravitating body, 



This is a priori obvious, because in this coordinate system, coordinate lengths don't depend on 
z. The Fermi coordinate system may be immediately derived from the Schwarzchild coordinates, 
far from the gravitating body (so that GM/r <C 1), by a simple coordinate redefinition. 

3.3. Alternative derivation 

For an alternative derivation of the force on the Casimir apparatus, we start from the variational 
definition of the stress tensor (I28|) . and consider a general coordinate transformation, 




(55) 



ds 2 = -(l + 2 5 z)dt 2 + dr 2 . 
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so that 



(58) 



where 



5g»u = Sx dxg^v + g 
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(59) 



For a rigid translation, 5x x is a constant, so only the first term here is present, which gives 
the result found above in section [34"1 However, if we don't make this restriction, we obtain 
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(60) 



where the surface term results from integration by parts. Again, notice for 5x x constant, the 
first and third terms identically cancel. If we take the region Q to be all space the surface term 
vanishes (and therefore so does the first term for a rigid displacement). A simple calculation 
indeed shows that the surface term cancels. 



3-4- Rindler coordinates 

Relativistically, uniform acceleration is described by hyperbolic motion 

t = ^sinhr, z = ^coshr, (61) 

which corresponds to the metric 

di 2 -dz 2 -dxi = £ 2 dr 2 - df - dx 2 ± . (62) 
The d'Alembertian operator has the following form in Rindler coordinates: 
d\ 2 ( d \ 2 „ 9 1 / d V 19 / J 



at) + UJ +vi = -e{a-r) + m\ e W + " i - m 

In this subsection, we will consider scalar fields and so-called semitransparent plates, that is, 
ones described by a <5-function potential 

V = XS(^-Q, (64) 

where £ a is the position of the plate. For a single semitransparent plate at £1, the Green's 
function can be written as 

= / ^(0 e "" (r " r ' )eik ' (M ' )±5(? '^ ) ' (65) 



where the reduced Green's function satisfies 



,2 



19 / J 



which we recognize as just the semitransparent cylinder problem with m — * C = ~' 1UJ an d 
k — > h [48]. The reduced Green's function for single plate thus is 

1 + A£i/ c (fc4i)A c (A;£i) 

= i c (ku )k c fc£> - \ " i , e,e'>ei- 67& 



In the limit A —* oo we recover the well-known Green's function subject to Dirichlet boundary 
conditions. 

If we use the uniform asymptotic expansion (UAE), based on the limit 

£ — > oo, £1 — ► oo, £ — £1 finite , ( = C£i — ► oo, (" finite , (68) 
we recover the Green's function for a single plate in Minkowski space, 

e — #*]€—€'] X e -«(l€-€i|+l€'-&l) 
fl«KiO -* A+~2k ' (W > 

where k = y/k 2 + £ 2 , = i£. 

The canonical energy-momentum for a scalar field is given by 

Tpv = d^d v 4> + gia/ -=C, (70) 



-5 

where the Lagrange density includes the ^-function potential. Using the equations of motion, 
we find the energy density to be 

From this and similar expressions for the spatial components of T^ u we can obtain {T^ v ) by 
applying a differential operator to the Green's function, in view of the one-loop connection (|25|) . 
The force density is given by 

fx = —^d v {^T\) + Ir^dxg^, (72) 

which follows from (|60j) . or 

/ € = -i^(^)-{T°°. (73) 
When we integrate over all space to get the force, 

F = J (dx)V=ff/{, (74) 



the first term in ([73]) is a surface term which does not contribute, as noted in (|60|) . so the 
gravitational force on the system is 



F = I = - I p T oo, (75) 



which when multiplied by the gravitational acceleration g is the gravitational force/area on the 
Casimir energy. Using the expression (|7ip for the energy density, and rescaling we see that 

dCd 2 fc ;2 



F= I d£C / tWC 2 5(£, (C = CO- (76) 



This result is an immediate consequence of the general formula [58] 

E c = jJ(dr) J ^2 W 2 S(r,r), (77) 



in terms of the frequency transform of the Green's function, 

^ e - M t-t> )g{r yy (7g) 

Alternatively, we can start from the following formula for the force density for a single 
semitransparent plate, 

/^I^A^-d), (79) 

or, in terms of the Green's function, 

T =- X - 2 \ -7^%[^(6,a)]- (80) 
For example, the force on a single plate is given by 

T= ~ d *k j ^T^ + WcW)]- ( 81 ) 
Expanding this about some arbitrary point £q, with ( = (£o, and using the UAE, we get 

T = 1 r <W ( . 2) 

QGttWo 1 + y/Aa' 1 ' 

which is just the negative of the (divergent) quantum vacuum energy of a single plate. 
For two plates at £i, £2, for < £1 < £2, 

u - T j, ^iKl + A 2 ^2^i - - ^1/2) rr ... , 

51?)? J — 1 < K > ^ [66a) 

while for £,£'>£ 2 >£i, 

) - /<*> - AAJ? + wl + hMi&iMMi - h«i) KKi m 

and between the plates, £1 <£,£'< £2, 

, K , = /<;f> _ A 2 fa*-j(l + A lWl)jJ) _ A 1& /?(l + A 2&gl / l)gg) 

where 

A = 1 + AiCiifi/i + A 2 6^2^2 - \iMtii2hK 2 {hK 2 - I 2 K{), (84) 

and where we have used the abbreviations I a = I((fc£ a ), / = I^(/c£), // = I^(k^'), etc. 

In the weak acceleration limit, the Green's function reduces to exactly the expected result, 
for £ x <£,£'< £ 2 (a = £ 2 -£i) 
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4k 2 

(85) 



with 
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1 + 



Ai 
2k 



1 + 



2k 



2ko, 



A1A2 

4k 2 ' 



(86) 



The flat space limit also holds outside the plates. 

In general, we have two alternative forms for the force on the two-plate system: 



T =-(%+%)i J 



1 f d(d 2 k 
"(2^" 



In A, 



(87) 



which is equivalent to ()76|) . From either of these two formulas, we find the gravitational force 
on the Casimir energy to be in the £ — ► 00 limit 

1 f°° 
T = — I d/c k 2 In Ao 
4?H Jo 



where Aq = e 2Ka A, A given in (|86f) . From this we get the explicit force 
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y+\ia ~ y+Ma 
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dyy' 
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Aia ^ 1 A 2 a ^ x 



-(5 C + £di + £42), 



which is just the negative of the Casimir energy of the two semitransparent plates, including 
divergent parts associated with each plate. The divergent terms simply renormalize the masses 
mi per area of each plate: 



£totai = mi + m 2 + £di + £d2 + £c = Mi + M 2 + £ c , 



(90) 



Mi = mi + £di , and thus the gravitational force on the entire apparatus obeys the equivalence 
principle 

gF = -g(M 1 + M 2 + £ c ). (91) 



3.5. Centripetal Acceleration 

Consider finally a Casimir apparatus undergoing centripetal acceleration with angular velocity 
u> (air <C 1), as illustrated in figure EJ A detailed calculation [59] shows that the force on the 
Casimir apparatus is just as expected, including the divergent pieces which renormalize the 
masses per area A of the plates: 

F uj 2 f 

J = -^ d 3 xrt 00 (r) = -uj 2 rcM (mi + m 2 + £ dl + £ d2 +£ c ), (92) 

where rcM is the center of energy of the entire system. Nontrivial is the fact that this result is 
independent of the orientation of the plates. 



4. Multiple Scattering Technique 

In the past few years, there has been a tremendous resurgence of interest in the multiple- 
scattering technique. This is largely inspired by the desire to improve the comparison with 
experiment, such comparison heretofore having been restricted to the use of the "proximity force 
approximation," which describes the Casimir interaction between bodies with curved surfaces by 
pairwise interactions between small plane elements of the surfaces [60]. The multiple-scattering 




Figure 3. Casimir apparatus undergoing uniform circular motion with angular velocity uj. The 
figure on the left has the Casimir plates tangent to the circle of revolution, while on the right 
the plates make an angle a with respect to that tangent. 

method certainly dates back to at least the Krein formula [61], and was first used in connection 
with the Casimir effect by Renne [62]. A most famous variation were the multiple reflection 
calculations of Balian and Duplantier [63]. In the last few years there has been an explosion of 
publications on variations of this technique, as physicists realized it could be practically used to 
obtain numerical results for forces between bodies of arbitrary shapes, for example [64-76]. 

We have also contributed to this redevelopment, giving the general formulation for scalar fields 
[77], and applying it to bodies described by 5-function potentials [78], and, after generalizing to 
electrodynamics, to dielectric bodies [79]. The methodology has also been applied to non-contact 
gears [80,81]. 

The multiple scattering approach starts from the well-known formula for the vacuum energy 
or Casimir energy (for simplicity here we first restrict attention to a massless scalar field) (r is 
the "infinite" time that the configuration exists) 

E = — TrlnG -> — TrlnGGn 1 
2r 2r 

where G (Go) is the Green's function, 

(-d 2 + V)G = 1, -d 2 G = l, 

subject to appropriate boundary conditions at infinity. We will use the Feynman or causal 
Green's functions. 

Now we define the T-matrix [82] 

T = S -1 = V(1 + G V)~ 1 . (95) 

If the potential has two disjoint parts, V = V\ + V 2 it is easy to derive 

T={V X + V 2 ){\ - G Ti)(l - GoTiGoTa)" 1 ^ - G T 2 ), (96) 

where 

T i = V i (l + G V i )-\ i = l,2. (97) 

Thus, we can write the general expression for the interaction between the two bodies (potentials) 
in two equivalent forms: 

E 12 = --^Trln(l - G Q T X G T 2 ) (98a) 
It 

= __L T rln(l - V^V^), (986) 



(93) 
(94) 



where G { = (1 + G Vi) Go, * = 1, 2. 



^.i. Multipole expansion 

To proceed to apply this method to general bodies, we use an even older technique, the multipole 
expansion. Let's illustrate this with a 2 + 1 dimensional version, which allows us to describe 
cylinders with parallel axes. We seek an expansion of the free Green's function 

„i|u| |r— R— r'| I' Ah 

G (R + r' - r) = J |r _ R _ r/| = j ^e^-*') 5Q(r± - R ± - r' ± ). (99) 

Here, the reduced Green's function is (u> = i() 

f (d 2 k\ \ p- ik i ,R ip ik i , ( r i- r i) 

As long as the two potentials do not overlap, so that we have rj_ — Rj_ — ^ 0, we can write 
an expansion in terms of modified Bessel functions: 

9o(rx - Rj_ - r'J = £ I^Kr^I^K^e-^f^nR), ^ = % + (*■ (101) 

m,m' 

By Fourier transforming, and using the definition of the Bessel function 

i m J m (kr)= / 27r ^ e - im< V fcrcos ^, (102) 
Jo 27r 

we easily find 

,0 , Km _ 1 f dfcfc J m (kr)J m (kr') (-IT' 

Thus we can derive an expression for the interaction between two bodies, in terms of discrete 
matrices, 

(t ^^T = lbl d ^ lndet i 1 ~ ~9°Tig 0T f 2 ) , (104) 
where T denotes transpose, and the T matrix elements are given by 

drrdcp [ dr'r , d07 m (Kr)e- im ^I m /(Kr')e im '^'T(r,0; r J , (/>'). (105) 



We will apply this 2+1 dimensional formalism to the interaction between two cylinders, as 
illustrated in figured Also illustrated there is the idea behind the proximity force approximation, 
which applies when the separation between the bodies is much smaller than the radii of curvature 
of the bodies. 

Consider two parallel semitransparent cylinders, of radii a and b, respectively, lying outside 
each other, described by the potentials 

Vi = X 1 d(r-a), V 2 = X 2 S(r' -b), (106) 

with the separation between the centers R satisfying R > a + b. It is easy to work out the 
scattering matrix in this situation, 

(*i)mm' = 2vrAia5 mm /-— — — — — — -. (107) 

1 + Xial m {Ka)K m (Ka) 



Figure 4. Geometry of two cylinders (or two spheres) with radii a and b, respectively, and 
distances between their centers of R > a + b. The proximity force approximation applies when 
the distance of closest approach d = R — a — b <C a,b. The approximation consists in assuming 
that the interaction is dominated by the interaction of adjacent surface elements, as shown. 



Thus the Casimir energy per unit length is 

1 f 00 

<B = — d/tKtrln(l — A), (108) 
47r Jo 

where A = B(a)B(b), in terms of the matrices 

BmnAa) = K m+m ,( K R) X \ aI }' { \^ . y (109) 

1 + \iaI m i(Ka)K m i{Ka) 

The strong coupling limit, the case of Dirichlet cylinders, has been considered in [75], for 
example. In contrast, here we will consider the situation of weak coupling, X±a <C 1, A26 <C 1, 
where the formula for the interaction energy between two cylinders becomes 

A A nh °° r°° 
{B = ~i^k- £ J te* K l + rn>Wl{™/R)ll>WR). (110) 



m,m'=— 00 



One merely exploits the small argument expansion for the modified Bessel functions I m (xa/R) 
and I m i(xb/R): 

40*0 = (2) E^M-n^) > ( 1U ) 
n=0 

where the coefficients Z mn are 

2 2(m+n) Y{m + n+^) 

m ' n ~ y^ n\ (2m + n)\ T(m + n + 1) ' ^ ' 

In this case we get an amazingly simple result 

n=0 

where fi = b/a, and where by inspection we identify the binomial coefficients 

= E (I) ^ fc - ( 114 ) 
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Figure 5. Plotted is the ratio of the exact interaction energy of two equal radii weakly-coupled 
cylinders to the proximity force approximation. 
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Figure 6. Plotted is the ratio of the exact interaction energy of two weakly-coupled cylinders 
with ratio of radii b/a = 99 to the proximity force approximation. 



Remarkably, it is possible to perform the sums, so we obtain the following closed form for the 
interaction between two weakly-coupled cylinders: 



A10A26 
8ttR 2 



a + b 
R 



R 



-1/2 



(115) 



We note that in the limit R — a — b = d ^ 0, d being the distance between the closest points 
on the two cylinders, we recover the proximity force theorem in this case 



U(d) 



The rate of approach is given by 



AiA 2 2ab 1 
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d <C a, b. 



(116) 
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R 2 -aR + a 2 d 
Aa{R-a) R' 



(117) 



The comparison between the exact result, and the proximity force approximation (|116|) is 
given in figure for equal radii, and in figure O for unequal radii. Evidently, the proximity force 
approximation fails badly as the cylinders are pulled apart. 
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Figure 7. Interaction energy between a weakly coupled (^-function cylinder and a Dirichlet 
plane. Plotted is the ratio between the energy and the PFA approximation (|12ip . 



We can similarly consider the interaction between a cylinder and a plane. By the method of 
images, we can find the interaction between a semitransparent cylinder and a Dirichlet plane to 
be 

1 f 00 

£ = — Kdnti ln(l - B(a)), (118) 
47r Jo 

where B{a) is given in (|109p . In the strong-coupling limit this result agrees with that given by 
Bordag and Nikolaev [72], because 

trB s = tvB s , B mm , = — ^ — -K m+ml (KR)I m ,{na). (119) 

K m {Ka) 

In exactly the way given above, we can obtain a closed-form result for the interaction energy 
between a Dirichlet plane and a weakly-coupled cylinder of radius a separated by a distance 
R/2. The result is again quite simple: 



Xa 



4vri? 2 



2a x 2 



-3/2 

(120) 



In the limit as d — > 0, this agrees with the PFA for this situation: 

Note again that this form is ambiguous: the proximity force theorem is equally well satisfied 
if we replace a by -R/2, for example, in U(d). The exact result is compared with the PFA in 
figure [7J 

4-2. 3- dimensional formalism 

The three-dimensional formalism is very similar. In this case, the free Green's function has the 
representation 

G (R + r / -r)= £ jKi|ClOip(i|CK)^(f)^(^Wm'(R)- (122) 

lm,l'm' 

The reduced Green's function can be written in the form 

<W™'( R ) = ( 4 -) 1 J ( 2vr)3^ + C 2 jKi|C|r)jHi|C|r') ?m( ^' m ' ( ^ ( } 



Now we use the plane-wave expansion 

e ik ' R = 47T Y, i'"il»(fcfl)^«m"(R-)^m"(k)» ( 124 ) 
l"m" 

so now we encounter something new, an integral over three spherical harmonics, 

J dkY^ m (k)y; m ,(k)^r, m „(k) = Q m ,iw,iv, (125) 

where 

lm,l m,l rn \ ) y ^ {0 J m' m" J K ' 

The three-j symbols (Wigner coefficients) here vanish unless I + I' + /" is even. This fact is 
crucial, since because of it we can write ji"(kR) in terms of Hankel functions of the first and 
second kind, using the reflection property of the latter, 

h$(kR) = {-lfh$(-kR), (127) 

and then extending the k integral over the entire real axis to a contour integral closed in the 
upper half plane. This leads to 



9h 



L.i'm'W = W '~V^ E Qm,Vrn',V'rn''K l/l+1/2 (\C\R)Y V/m/l (R). (128) 



l"m" 

For the case of two semitransparent spheres that are totally outside each other, 

V l (r) = X l 5(r-a), V 2 (r') = \ 2 5{r' - b), (129) 

in terms of spherical coordinates centered on each sphere, it is again very easy to calculate the 
scattering matrices, 

and then the harmonic transform is very similar to that seen for the cylinder, (A; = i|C|) 
(h)im,l'm> = J (dr)(dr')iK^)^m(f)jH^ , )^'m'(f / )ri(r,r / ) 

X A f i\l XlQ7T I l+i/ 2 (\C\ a ) 

" " mm ' [ ' 2|C| I + XiaK l+1/2 (\C\a)I l+1/2 (\C\aY ^ j 

Let us suppose that the two spheres lie along the z-axis, that is, R = Rz. Then we can simplify 
the expression for the energy somewhat by using Y[ m (Q = 0) = <5 m o\/ (2/ + l)/47r. The formula 
for the energy of interaction becomes 



E = — 

2vr 



^°°dCtrln(l-A), (132) 



in terms of the matrix 

Alm,l'm> = ^m,ra' '/ J Bii" m (a)Bi"lf m (b) (133) 
I" 

with 



B Wn {a) = -^glr'+V(2i + mi' + 1) E(2/" + 1) 



ll'l"\fl V l"\ ^''+i/ 2 (C^)Aia/ / 1 +1/2 (Ca 



00 O^m-m 0; I + \iaI v+1/2 {Ca)K v+l/2 {Ca) 



(134) 



Note that the phase always cancels in the trace. 

Again, the strong coupling limit was considered in [75]. For weak coupling, a major 
simplification results because of the orthogonality property, 

£ (m -m ) (m -m ) = 5v ' l '"WTl' 1 ~ (135) 

m=-l V 7 V 7 

and then we find for the interaction energy 

^_^^| E(21 + 1)(21 , + 1)(2i „ + 1)( uw : ^ 

• yu ll'l" v 7 

xKf fl+1/2 (x)lf +1/2 (xa/R)lf >+1/2 (xb/R). (136) 

As with the cylinders, we expand the modified Bessel functions of the first kind in power 
series in a/R,b/R < 1. This expansion yields the infinite series 
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\ ia \ 2 bab^ 1 ^ /a\2(n-m)/b\ 2m 



/. >2 ,., „ 

n=U m=0 



where by inspection of the first several D n ^ m coefficients we can identify them as 

1 / 2n + 2 \ . . 

D -™ = 2{2 m + 1 } (138) 

and now we can immediately sum the expression for the Casimir interaction energy to give the 
closed form 



^^W^l. (139) 



< a+b \ 2 
\ R J 

Again, when d = R — a — b <C a,b, the proximity force theorem is reproduced: 



U(d) ~ ln(d/R), d « a, 6. (140) 

167m 

However, as figure [8] and figure [9] demonstrate, the approach is not very smooth, even for equal- 
sized spheres. The ratio of the energy to the PFA is (b/a = /x) 

E Mil + u) 2 /2u] , , . 

U- 1+ \J/R • "«"'"■ (141) 
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Figure 8. Plotted is the ratio of the exact interaction energy of two weakly-coupled spheres of 
equal radii to the proximity force approximation. 
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Figure 9. Plotted is the ratio of the exact interaction energy of two weakly-coupled spheres, 
with ratio of radii b/a = 49, to the proximity force approximation. 




In just the way indicated above, we can obtain a closed-form result for the interaction energy 
between a weakly-coupled sphere and a Dirichlet plane. Using the simplification that 

£(-ir(^Jol(ooo)=*»> < 142 > 

m=-l v / \ / 

we can write the interaction energy in the form 

E = -^kl ^E^< 2, + 1)/ ^w^(4)' (143) 

Jl) i=o v 

Then in terms of R/2 as the distance between the center of the sphere and the plane, the exact 
interaction energy is 

A /a\2 1 



e = -yAr) i-(2./«f < 144 » 

which as a — * R/2 reproduces the proximity force limit, contained in the (ambiguously defined) 
PFA formula 

""S3- (145) 

The ratio of the exact interaction to the PFA is plotted in figure [101 
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Figure 10. Plotted is the ratio of the exact interaction energy of a weakly-coupled sphere above 
a Dirichlet plane to the PFA. 

4-3. Exact Results — Weak Coupling 

In weak coupling it is possible to derive the exact (scalar) interaction between two potentials in 
general forms, 



r — r 



3 D: £ = -^/(dr)(drO^|^. (1466) 

The results found above in sections 14.11 and 14.21 can be directly derived from these pairwise 
summation formulae. For exact results for (^-function plates, see the contribution to this 
Proceedings by Jef Wagner. 

It is straightforward to adapt the multiple scattering formalism to the far more relevant 
situation of electromagnetism. For weak coupling, this reduces to the (retarded dispersion) van 
der Waals (vdW) potential between polarizable molecule given by 

23 aiot 2 e-l 

This allows us to consider in the same vein (electromagnetic) interactions between distinct dilute 
dielectric bodies of arbitrary shape. (Of course, this potential represents the retarded limit of 
the van der Waals force. The non-retarded limit of the van der Waals force between bodies has 
long been studied, for example between dielectric cylinders [83,84], and between a sphere and a 
plane [85].) 

This vdW potential may be directly derived from the multiple scattering formula 

W = -L Tr l n JL~-l TrVl r V 2 T , (148) 
l l o l 

where where have omitted self-action terms, and 

e -|CI|r-r'| 

r = VxVxl— - + (149) 



The quantity V = £ — 1 is the potential representing the dielectric body. 

We first apply this approach to compute the force between a dilute dielectric slab and an 
infinite dilute dielectric plate, as shown in figure fTTl If the cross sectional area of the finite slab 
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Figure 11. A plate of cross section A, a distance a above an infinite plate. Both plates are of 
semi-infinite extent in the direction perpendicular to both plates. 



is A, the force between the slabs is 

F 23 1 



A 640vr 2 a 4 



(ei - l)(e 2 - 1), (150) 



which is the Lifshitz formula for infinite (dilute) slabs. Note that there is no correction due to 
the finite area of the upper slab. 

In the same way we can immediately compute the force between a dilute dielectric sphere 
of radius a a distance Z away from an infinite slab, as illustrated in figure [T2"l The interaction 
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Figure 12. A dilute dielectric sphere, of radius a, the center of which is a distance Z above an 
infinite dilute dielectric plate. 

energy of this system is 

23 . w ,4vra 3 /3 1 
£= -640^' e '- 1 >' e2 -" Z' <151) 

which agrees with the PFA in the short separation limit, Z — a = 5 <C a, where the force is 

F PFA = 2ira£ l{ (6) = -gj^i - l)(e 2 - 1)^, (152) 

with an exact correction, intermediate between that for scalar averaged Dirichlet and Neumann 
[71,72] and electromagnetic perfectly-conducting boundaries [73]. 
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Figure 13. Rectangular solid (with dielectric constant £2) of side a, 6, and c (perpendicular to 
the plane, not shown) a distance Z above an infinite plate (with dielectric constant e\) extending 
to z = —00. The side b makes an angle 9 with respect to the plate. 
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Figure 14. The stable equilibrium angle of a slab above an infinite plate for given b/a ratios 
0.95, 0.9, and 0.7, respectively given by solid, dashed, and dot-dashed lines. 




We can also consider torques between bodies. For example, consider a rectangular dilute 
dielectric solid, with sides a, b, and c, the center of mass of which is a distance Z above a dilute 
dielectric plate, but which is tilted through an angle 9 relative to the plate, as shown in figure [T3l 
Generically, the shorter side wants to align with the plate, which is obvious geometrically, since 
that (for fixed center of mass position) minimizes the energy. However, if the slab has square 
cross section, the equilibrium position occurs when a corner is closest to the plate, also obvious 
geometrically. But if the two sides are close enough in length, a nontrivial equilibrium position 
between these extremes can occur, as well illustrated in figure [TH For large enough separation, 
the shorter side wants to face the plate, but for 



„ _ a / 2a 2 + 5b 2 + V9a 4 + 20a 2 6 2 + 206 4 , , lco , 

Z<Z = -^ , a>b, (153) 

the equilibrium angle increases, until finally at Z = D = \/« 2 + b 2 /2 the slab touches the plate 
at an angle 9 = arctan6/a, that is, the center of mass is just above the point of contact, about 
which point there is no torque. 

Finally, once again we can get exact results for the vdW interaction between parallel cylinders 



(radii a, b, center separation R): 
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and between spheres: 
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(155) 



This last expression, which is rather ugly, may be verified to yield the proximity force theorem: 

23 a(R-a) 



E^U 



640vr R5 2 



5 = R — a — b <C a,b. 



(156) 



It also, in the limit b — > oo, R — > oo with R — b = Z held fixed, reduces to the result (|15ip for 
the interaction of a sphere with an infinite plate. 

The same methodology may be used to derive the force between corrugated planes and 
corrugated concentric cylinders, which calculations may have relevance to nanomechanics. See 
the contribution to this Proceedings by Cavero-Pelaez, and [80,81]. 



5. Comments and Prognosis 

In this talk I have described some recent exciting progress in the practical development of 
quantum vacuum energy. I have concentrated on the work of the Oklahoma group, without 
meaning to denigrate the work of many others of equal or greater importance. I will summarize 
the status of the topics I discussed briefly here: 

• The temperature controversy can only be resolved by dedicated experiments, designed to 
measure Casimir forces at different temperatures between metal surfaces 

• The semiconductor anomaly needs to be examined more closely both theoretically and 
exp er iment ally. 

• Casimir energies, including their divergent contributions, exhibit the inertial and 
gravitational masses expected: m = E/c 2 . These gravitational interactions need to be 
examined for more complicated geometries, such as for the Casimir energy of a perfectly 
conducting spherical shell, with the hope of understanding the physical significance of the 
divergent contributions. 

• The multiple-scattering methods recently developed are in fact not particularly novel, and 
illustrate the ability of physicists to continually rediscover old methods. What is new is the 
recognition that one can evaluate continuum determinants (or infinitely dimensional discrete 
ones) accurately numerically, and in some cases even exactly in closed form. This is making 
it possible to compute Casimir forces and torques for geometries previously inaccessible. 

• It is indeed remarkable, if perhaps not surprising in retrospect, to see that closed form 
expressions can be obtained for the interaction between spheres, between parallel cylinders, 
and for other geometries in weak coupling. These results demonstrate most conclusively the 
unreliability of the proximity force approximation (of course, the proximity force theorem 
holds true). 



• This methodology has been used to obtain new results for non-contact gears: See the talk 
by Ines Cavero-Pelaez. 

• Jef Wagner has presented new results for edge effects for semitransparent potentials. 

• Undoubtedly these developments will lead to improved conceptual understanding, and to 
better comparison with experiment. 
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